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Abstract 

The topological susceptibility is one of the few physical quantities that directly measure 
the properties of the QCD vacuum. Chiral perturbation theory predicts that in the small 
quark mass limit the topological susceptibility depends quadratically on the pion mass, 
approaching zero in the chiral limit. Lattice calculations have difficulty reproducing this 
behavior. In this paper we study the topological susceptibility on dynamical staggered 
fermion configurations. Our results indicate that the lattice spacing has to be small, 
around a ~ O.lfm for thin link staggered fermion actions to show the expected chiral 
behavior. Our preliminary result indicates that fat link fermions, on the other hand, 
reproduce the theoretical expectations even on lattices with a ^ 0.17fm. We argue that 
this is due to the improved fiavor symmetry of fat link fermionic actions. 
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I. INTRODUCTION 



Instantons play an important role in the QCD vacuum from the breaking of the ax- 
ial Ua{^) symmetry to chiral symmetry breaking and the low energy hadron spectrum||l|. 
Lattice studies support many of the theoretical predictions^. For example, the Witten- 
Veneziano formula relates the pure gauge topological susceptibility to the masses of the 
?7, 7]' and K mesons predicting xH'^ = ISOMeV. The topological susceptibility on pure gauge 
lattices has been measured by several groups using different methods obtaining consistent 
results xj/' = 205(5)MeVi, |, §. 

Light fermions suppress the topological susceptibility. Chiral perturbation theory 
predicts [|] 

x = ^s + oK) = ^ + oK), (1) 

where ruq is the quark mass, is the number of fermion flavors and E is the chiral conden- 
sate per fermionic flavor. In this normalization = 92MeV. Several recent lattice studies 
measured x oii dynamical conflgurations. Calculations with Wilson and clover fermions cover 
a fair range of lattice spacing a ^ 0.08 — 0.20fm and pion masses m^ro ~ 1.3 — 2.5 (ro is the 
Sommer parameter 0) |||, H, |10|, |ll], The results appear controversial. UKQCD uses 
clover fermionslH, ||, ITof. Their data on lattices with a ~ O.lfm is basically consistent with 
eqn.(ll). SESAM/Tx-^^ uses unimproved Wilson fermions. Their topological susceptibility 
at similar lattice spacing does not decrease with the pion mass though the statistical errors 
are too large to claim inconsistency with theoretical expectations [0]. CP-PACS published 
data using clover fermions and improved gauge action at lattice spacing a ~ 0.17fm|[T2|. 
Their conclusion is the same as SESAM/T^L. The situation with staggered fermions is not 
much better. Only the Pisa group measured the topological susceptibility with two and four 
flavors of staggered fermions at lattice spacing a ~ 0.09 — 0.17fm[|l^, |14|. They do not see 
the reduction of the topological susceptibility with decreasing quark mass either. 

Can we understand what is going on with the dynamical simulations? Eqn.([l|) is valid 
only in case oiuf light chiral fermions creating n^f — ^ light pions. Both Wilson and staggered 
fermions violate this assumption. Wilson fermions break chiral symmetry explicitly. The 
addition of the clover term reduces the symmetry breaking and improves chiral behavior. 
Staggered fermions have only a residual U{1) chiral symmetry and only one true Goldstone 
boson, the other pion-like states can be heavy. Fat link fermions considerably improve flavor 
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symmetry and consequently chiral behavior. Both fermionic formulations should show the 
expected chiral behavior in the continuum limit, clover and fat link fermions sooner than the 
unimproved ones, but it is not clear when scaling in the topological susceptibility actually 



sets in. In a recent paper |^] it was suggested that the reason the topological susceptibility 
from available lattice data does not follow the expected theoretical behavior is a combination 
of three effects: too large lattice spacing, too small volumes and too small Leutwyler-Smilga 
1^ parameter. While all three conditions could indeed be important, we feel that the non- 
chiral behavior of the fermionic actions is the main cause of the problem. Chiral symmetry 
violation of the fermionic action is a scaling violation effect and is covered by the first 



condition of ||T5|. However, it is not a universal quantity, it can strongly depend on the 
fermionic action. 

Why are the topological properties of the vacuum important? Phenomenological in- 
stanton models predict that the low energy hadron spectrum is governed by the near-zero 
eigenvalue modes of the Dirac operator, which, in turn, are related to instantons. If the 
fermion-instanton interaction is different from the continuum one either because of chiral 
symmetry breaking or other lattice artifacts, the Dirac spectrum and consequently the low 
energy hadron spectrum can also be different. Recently, using chiral symmetric overlap 
fermions on the lattice, we showed that with light quarks the first few modes of the Dirac 
operator saturate the low lying hadron propagators on quenched a ~ 0.12 lattices, just as 



the phenomenological models predict ||TB|, O. In a subsequent publication contradictory con- 



clusion was reached using Wilson fermions ||l8l. For us that implies that Wilson fermions, at 
least at large to moderate lattice spacings, interact differently with instantons than chirally 
symmetric fermions. 

In this paper we investigate the topological susceptibility on both Uf = 2 and nj = i 
staggered dynamical fermion configurations. The Uf = 2 configurations are 16^32 thin link 
staggered fermion lattices created by the Columbia and MILC collaborations and down- 
loaded from the NERSC archive ||19|. The Uf = A lattices are smaller, 8^24 configurations 



used in the study of fat link fermions in Ref. |^. Two of the n/ = 4 sets were generated 
with thin link staggered fermions and one with N=3 level APE blocked fat link fermions. 
The latter action has about an order of magnitude better flavor symmetry than thin link 
actions at similar parameter values. The thin link staggered fermion results with a < O.lfm 
are more or less consistent with eqn.(fD. The thin link a ^ 0.1 7fm data shows clear devi- 
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ation, the topological susceptibility is consistent with the quenched value, independent of 
the pion mass. On the other hand, the fat link data at the same coarse lattice spacing is 
in perfect agreement with eqn.(jl|), suggesting that improved flavor/chiral symmetry indeed 
has a strong effect on the topology. 

To determine the topological charge of the configurations we used a topological charge 
density operator constructed from hypercubic blocked (HYP) fat links||21|. Hypercubic 
blocking was introduced in a recent paper as an alternative to repeated APE smearing ||22|. 
Hypercubic blocking mixes links only from the hypercubes that attach to a given link, thus 
the fat link is very compact yet the configuration is as smooth as after three levels of APE 
blocking. To avoid the distortion effect of extended operators, we consider only one to three 
levels of hypercubic blocked operators (HYPl, HYP2 and HYPS). We have calculated the 
additive and multiplicative renormalization factors of the topological susceptibility for these 
operators. After two to three levels of HYP blocking the renormalization factors turn out 
to be consistent with their tree level values and, after correction, all three HYP topological 
susceptibility measurements are consistent. 

In Sect. 2. we describe the hypercubic topological operator and illustrate the measure- 
ment of the renormalization factors on pure gauge Wilson plaquette configurations. Sect. 3 
contains our results for two- and four-fiavor staggered fermions. Sect. 4 is a short discussion 
on fiavor symmetry and the summary of our results. 

II. HYPERCUBIC TOPOLOGICAL CHARGE OPERATOR 

Most large scale simulations use pure gauge observables to measure the topological sus- 
ceptibility with some discretized version of the continuum FF as the lattice charge density 
operator qiix). The relation between the continuum and lattice topological susceptibilities 
contains both a multiplicative and an additive renormalization factor. The lattice charge 
density operator is related to the continuum one through a multiplicative renormalization 
factor 



In addition the correlator of two topological density operators have an additive correction 
term as well due to the mixing of qiix) with other lattice operators 



qiix) = Za^q{x). 



(2) 



QL{x)qL{0) = Z^a^g(x)g(0) +m{x), 



(3) 
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thus the topological susceptibility on the lattice is 

XL = J q{x)q{0)d'x = Z^a\ + M, (4) 
M = [ m{x)ct^x. (5) 



The renormalization factors Z and M depend on the lattice parameters (3 and m^. The renor- 
malization constants in principle can be determined non-perturbatively using the heating 
method proposed by the Pisa group||23|. 

Due to lattice artifacts, mainly dislocations, on typical lattices Z = 0.25 and M/xl — 1 
for thin link topological density operators which makes the direct determination of x almost 
impossible. Local smearing or cooling removes most of these lattice artifacts moving Z 
towards 1 and M to 0. Repeated smoothing gives Z = 1 and M/x -C 1 leading to Xl — ci^X- 
Since the non-perturbative calculation of Z and M introduces statistical and systematical 
errors, a topological density operator where the renormalization constants are small or can 
be neglected is desirable. However, repeated smearing and cooling methods have their 
drawbacks as well. While removing vacuum fluctuations and lattice artifacts, both methods 
remove topological objects, mainly small instantons and close-by pairs. In addition the size 
of the remaining objects change as well, as can be demonstrated by monitoring individual 
instantons during the smearing process^. 

In this paper we construct the topological charge density operator qiix) using hypercubic 
blocked links and the improved thin link charge operator of Refs. |^ and Hypercubic 



blocking was introduced in Ref. ||21| as a localized alternative to repeated APE blocking. 
HYP fat links mix original links from the hypercubes that are attached to the fat link only, 
yet they create configurations that are as smooth as the ones obtained after three levels of 
APE blocking. 

Topological charge measurement are difficult because of the presence of dislocations: short 
distance vacuum fiuctuations that can be mistaken for small instantons. Smearing the links 
in the topological density operator reduces this problem in part by sharpening the transition 
between the topologically non-trivial Q = J d'^xq{x) = 1 sector and the trivial Q = sector. 
Figure [l| shows the topological charge of a smooth instanton measured with the thin link, 
HYPl, HYP2 and HYPS fat link operators as the function of the instanton radius . The 
charge measured with the HYP operators rises sharply at instanton radius p/a ~ 1. 
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Figure 1: The topological charge of a smooth instanton as the function of its radius pja. The 
charge is measured with thin link operator (crosses) , HYPl (diamonds) , HYP2 (bursts) and 
HYPS (octagons) fat link operators. 
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Figure 2: The topological charge on 500 heated 8'^ Q = configurations measured with a) HYP2 
topological operator, b) HYPl topological operator. The configurations were heated with 10-50 
heat bath steps using /3 = 6.0 Wilson pure gauge action. 



Even though the HYP operators remove most dislocations, they still can have non-trivial 
renormalization factors. We have measured the renormalization factors Z and M following 



the heating method of Ref. p^]. To measure M one has to heat a trivial Q = configu- 
ration. The configuration thermalizes fast at short distances, and since the origin of M is 
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Z 


MxIQS 


XL X 10^ 


a^X X 10^ 


4 


APEl 


0.795(4) 


0.95(7) 


5.4(5) 


7.0(7) 


0.058(6) 


HYPl 


0.935(4) 


0.11(3) 


6.3(4) 


7.1(5) 


0.059(6) 


HYP2 


1.000(4) 


0.006(4) 


7.0(4) 


7.0(4) 


0.058(6) 



Table I: Results for the renormalization factors Z and M and the lattice and continuum topological 
susceptibilities for (3 = 6.0 pure gauge Wilson action. The topological susceptibility was measured 
on 220 16'^32 configurations downloaded from the NERSC archivep5||. 

local contact terms, short distance thermalization is sufficient to measure M. Before heating 
creates non-trivial topological objects the thermalization must be terminated and restarted 
with a different random seed. It is important to make sure the measurement of M is done 
on Q = configurations only. Fortunately with the HYPS and frequently with the HYP2 
operators it is easy to separate the trivial configurations from the occasional Q 7^ ones. 
Figure ^a shows the topological charge on 500 heated (5 = configurations measured with 
HYP2 operator. Ordered 8"^ configurations were heated with 10-50 heat bath steps using 
13 = 6.0 pure gauge Wilson action. Most configurations have very small topological charge, a 
few has ~ 1 and only 2-3 configurations have topological charge whose interpretation is 
not clear. I chose, based on Figure |l| but somewhat arbitrarily, a cut Qmax = 0.3 to separate 
the Q = and Q 7^ configurations. The final results are fairly insensitive to the precise 
choice of this cut. The topological charge measured with the HYPl operator fiuctuates more 
as is illustrated in Figure |]/b. Again, a configuration is accepted as Q = if its topological 
charge is \Ql\ < 0.3 as measured with the HYP2 operator. The multiplicative renormaliza- 
tion constant Z is calculated similarly on Q = 1 configurations that contain a single smooth 
instanton of size p/a = 3.0,2.5 or 2.0. Z can be different on different instanton size back- 
grounds, though in none of the measurements we performed was the difference more than a 
few per cent. We chose Z on the background where the smooth instanton's size was closest 
to the expected average instanton size of the Monte Carlo configurations, p/a ^ 0.3fm/a. 

Results for /3 = 6.0 pure gauge Wilson action are summarized in Table |I| where Z, M, 
Xl, X(^^ = {Xl — M)/Z'^ and x'^q with Sommer parameter tq = 5.37(1) are given both 
for the HYPl and HYP2 operator, and, for comparison, one level APE smeared operator 
with parameter a = 0.75 as well. The topological susceptibility was measured on 220 16^32 
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Z 


MxlO'' 


XL X 10^ 


a^X X 10'' 


4 


HYPl 


0.85(4) 


0.6(1) 


5.1(9) 


6.2(14) 


0.048(12) 


HYP2 


0.95(1) 


0.19(4) 


6.0(9) 


6.4(12) 


0.049(11) 


HYP3 


0.98(1) 


0.02(1) 


6.2(9) 


6.4(10) 


0.049(9) 



Table II: Results for the renormalization factors Z and M and the lattice and continuum topolog- 
ical susceptibilities for the (i = 5.7 pure gauge Wilson action. The topological susceptibility was 
measured on 350 8^24 configurations. 

/? = 6.0 configurations from the NERSC archive ||25i. The results for a^x c^re consistent for 
all three operators, in physical units x^^^ = 196(5)MeV. The renormalization factors, on 
the other hand, are quite different. The background term M is about 17% of xl for the 
APEl operator, 2% for the HYPl operator and less than 0.1% for the HYP2 operator. The 
renormalization factors for the HYP2 operator are negligible, for all practical purposes we 
can use Z = 1 and M = 0. 

Results are similar at (3 = 5.7 though there one needs three levels of HYP blocking to 
reduce M to and Z to 1 as is illustrated in Table |I|. Here the topological susceptibility 
was measured on 350 8^24 lattices and the renormalization constants are obtained from 500 
heated 8^^ lattices. The predictions for the continuum topological susceptibility are consistent 
from all three operators, x^^^ = 188(9)MeV. Within statistical errors it is also consistent, 
though about 15% lower, then the value obtained on the /3 = 6.0 data set. Whether the 
difference is only statistical fiuctuation or the a ^ 0.17fm lattice spacing at /? = 5.7 is too 
coarse to support all the topological objects on the lattice cannot be decided with the present 
statistics. 

III. THE TOPOLOGICAL SUSCEPTIBILITY ON DYNAMICAL CONFIGURA- 
TIONS 

A. Two flavors of staggered fermions 

We have analyzed several two-fiavor configuration sets using HYP1-HYP3 operators. 
The renormalization factors have to be calculated independently at every parameter value. 
The results for both Z and M were very similar to the quenched values at similar lattice 
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^conf 


ro/a 


a[fm] 




a^X X 10^ 


(rom^)2 


4 


13 = 5.7,amq = 0.01 


83 


6.29(8) 


0.08 


0.252(2) 


2.6(3) 


2.5(1) 


0.041(6) 


f3 = 5.7, om, = 0.015 


46 


6.02(8) 


0.08 


0.293(2) 


3.5(6) 


3.1(1) 


0.046(10) 


f3 = 5.7, amg = 0.025 


33 


5.8(1) 


0.09 


0.388(1) 


5.8(9) 


5.1(2) 


0.065(15) 


(3 = 5.415, amg = 0.025 


201 


2.96(3) 


0.17 


0.4454(2) 


84(8) 


1.74(1) 


0.064(9) 



Table III: Results for the topological susceptibility on = 2 configurations. All lattices are 16^32 
standard thin link staggered fermion configurations. 

spacing. We found that on configurations with lattice spacing a ~ O.lfm the renormalization 
constants of the HYP2 operator could be neglected, while configurations with a ^ 0.17fm the 
HYP2 operator had a few percent correction from the renormalization constants, the HYPS 



operator had none. Table III collects our results. All four data sets are from the NERSC 



archive, the first three were generated by the Columbia group, the last one by MILC|2^, |27 . 
All lattices are 16^32 and use standard thin link staggered fermions with plaquette Wilson 
gauge action, am^ was measured in the original studies, ro/a for the MILC lattice is from 
Ref. [|^. For the Columbia lattices we have measured rg/a using HYP blocked Wilson 
loops. The HYP potential has greatly reduced statistical errors making it possible to obtain 
reliable values even from 33-83 configurationsp|]. The lattice spacing in the fourth column 
was obtained using ro = 0.5fm and is listed for future reference. Since the different HYP 
topological charge operators give consistent results, only the continuum value a'^x is listed. 
The topological susceptibility for the first data set has been measured in Ref. [^] using 10- 
40 APE smeared operators. The values for a'^x with APE smeared operators are consistent 
with the present result. In Figure H/a x'^o i^ plotted as the function of (m^ro)^. The filled 
octagons correspond to the Columbia data sets with a < O.lfm. The filled diamond is the 
MILC lattice data with a ~ 0.17fm. The filled square at the lowest m^ro value is also 
from MILC. It is on a 24^64, /? = 5.6, amg = 0.01 configuration set and was measured by 



C. DeTar using 10-20 level APE blocked topological operators]^. The lattice spacing on 
these configurations is a ~ O.lfm, ro/a = 4.99. The dashed line in Figure ^ is the leading 
order theoretical prediction from eqn.(|) using fj, = 92MeV. The three lowest m.„rQ data 
points on the finer lattices with a < O.lfm are consistent, though a bit higher then the 
leading order theoretical curve. However, the data point with a ~ 0.17fm is very different. 
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Figure 3: x^t the function of {m.,^rQ)^ with nj = 2 staggered fermions. The dashed hne is the 
leading order theoretical prediction from eqn.([l]) using f.,^ = 92MeV. The dotted lines on the right 
indicate the quenched value of x'^o- ^) Staggered fermions. Filled octagons: Columbia data sets; 
filled diamond: MILC data set; filled square: MILC||3^. b) Same as a) with Wilson/clover data 
added. Octagons: UKQCD; diamonds: SESAM/TxL. 



The topological susceptibility on those lattices is consistent with the quenched result even 
though the pion mass is fairly small. Since on quenched configurations the topological 
susceptibility can be measured successfully even on coarse lattices, this discrepancy is not 
likely to be the consequence of the gauge system. Rather, it appears that the instantons do 
not feel the effect of light staggered fermions on coarse lattices. Flavor symmetry violation 
of staggered fermions can explain these findings. In Figure ^b published Wilson/clover 
data is added to the a < O.lfm staggered data. Octagons correspond to the clover fermion 
simulations of UKQCD, diamonds to the Wilson fermion simulations of SESAM/T^L- All 
Wilson/clover data has lattice spacing a < O.lfm. The data from UKQCD is consistent 
with the staggered results and with the leading order theoretical curve as well. However the 
topological susceptibility with unimproved Wilson fermions does not follow the expected 
behavior at small pion mass. It remains large, consistent with the quenched result, the same 
behavior we saw with staggered fermions at lattice spacing a ~ O.lTfm. It appears that 
Wilson fermions have a very different effect on instantons even on finer, a ~ O.lfm lattices 
than clover or even unimproved staggered fermions. This is likely the consequence of chiral 
symmetry breaking of the Wilson action. 
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ro/a 


a[fm] 


arriT^ 


a\ X 10'' 


(rom^)2 


4 


f3 = 5.2, am, = 0.06, thin 


188 


2.85(2) 


0.18 


0.661(1) 


8.5(9) 


3.5(1) 


0.056(8) 


P = 5.25, aniq = 0.06, thin 


189 


3.48(3) 


0.14 


0.664(1) 


4.8(5) 


5.3(1) 


0.070(10) 


f3 = 5.2, anig = 0.1,fat 


140 


2.97(3) 


0.17 


0.695(4) 


3.5(4) 


4.2(2) 


0.026(4) 



Table IV: Results for the topological susceptibility on n/ = 4 configurations. All lattices are 8^24. 
The first two sets were generated with standard thin link staggered fermions, the third one with 
N=3 APE blocked fat link fermions. 

B. Four flavors of staggered fermions 



Unfortunately we could not find any large n/ = 4 staggered data sets to use. In Ref. ||20 
where we proposed a dynamical fat link fermion update, we generated three 8^24 n/ = 4 
data sets to study flavor symmetry breaking of thin and fat link actions. The flrst two 
sets were generated with standard thin link staggered fermions, the third one with N=3 
APE blocked fat link fermions. All three sets have lattice spacing a ^ 0.17fm as shown in 
Table 0. The pion mass values in the table are also from Ref.|^ but we have re-analyzed 
the potential data using HYP Wilson loops and the ro/a values listed in Table |V| are more 



reliable and slightly different from the ones used in Ref. pO| The difference in the topological 
susceptibilities between the thin and fat link actions is striking. There is no Wilson/clover 
fermion data available at = 4 to compare the staggered result. In Figure || we plot x^t 
as the function of {mT,rQY for the three data sets from Table The dashed line is the 
leading order theoretical prediction from eqn.(|l|) using the experimental value /,r = 92MeV. 
The fat link action data shows the expected behavior while the thin link action topological 
susceptibility values are consistent with the quenched result, independent of the dynamical 
quark mass. For the thin link action this is the same behavior we saw with the nj = 2 data 
at similar lattice spacing. 



IV. DISCUSSION AND SUMMARY 



The interaction between light quarks and instantons changes the QCD vacuum substan- 
tially. One of the easiest way to get information about this effect is through the topological 
susceptibility that is expected to scale with the square of the pion mass in the small quark 
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Figure 4: xi^o as the function of (m,r'^o)^ with = 4 staggered fermions. The filled diamonds are 
from thin link staggered configurations, the filled octagon is from fat link dynamical configurations. 
The dashed line is the leading order theoretical prediction from eqn.(|l|) using the experimental 
value /tt = 92MeV. The dotted lines on the right indicate the quenched value of x^o- 

mass limit. 

Our results indicate that thin link staggered fermions at lattice spacing a > 0.1 7fm 
do not have the expected continuum effect on instantons, a lattice spacing a ~ O.lfm is 
needed to recover the proper chiral behavior. In contrast to that, fat link staggered fermions 
show the expected behavior even on coarser a ^ 0.17fm lattices. Can this be understood 



as the effect of flavor symmetry violation? In Refs. ||20|, |21| we studied flavor symmetry 
violation on quenched lattices both with the thin link, N=3 APE blocked and hypercubic 
blocked fat link fermions. The two fat link formulations have about the same level of 
flavor symmetry violations as measured by the relative mass splitting = (m^r — mG)/rnG 
between the Goldstone pion and the other pion like objects, m.^-. For thin link fermions 



at (m^ro) 



2.0 we found ^ 0.7 at a ^ 0.17fm, A^ ^ 0.2 at a ^ O.lfm for the 



lightest non-Goldstone pions. The corresponding values for fat link fermions are A^ ^ 0.09 
and A^ ^ 0.01, respectively. QCD with two(four) light flavors should have 3(15) light 
pions. Apparently when there is only one light pion and the other pion-like objects are 70% 
or more heavier than the Goldstone pion, the vacuum does not look like a two (four) -flavor 
QCD vacuum. The flavor symmetry breaking has to be reduced below 20% to get acceptable 
results. Fat link fermions can do that easily even on coarse lattices. 
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Available lattice data for Wilson/cover fermions can be understood similarly. The topo- 
logical susceptibility with clover fermions are consistent with the theoretical predictions at 
lattice spacing a ^ O.lfm but the Wilson fermion data at the same lattice spacing indicates 
that unimproved Wilson fermions, that have much larger chiral symmetry violations, do not 
have the expected continuum like interaction even at a ~ O.lfm. Results from CP-PACS 
at lattice spacing a ^ 0.17fm indicate that at that lattice spacing not even clover fermions 
can reproduce the continuum topological behavior. Fat link clover fermions improve chiral 
symmetry and could provide a better alternative to thin link clover fermions. It would be 
interesting to find a parameter similar to A^r that correlates with the level of chiral symmetry 
breaking for the Wilson like fermions and compare it with the topological susceptibility. 

One should be concerned about dynamical simulations where the topological susceptibility 
is not reproduced correctly, since that indicates that the vacuum at those simulations is 
more like the quenched vacuum rather than the expected dynamical one. Improving chiral 
symmetry can have a profound effect. This point is further underscored in a forthcoming 



publication about the finite temperature phase diagram obtained with fat link fermions|31 
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